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Imaging 

This invention relates to imaging and in particular to 
5 methods and systems for image enhancement. 

Imaging is carried out by transfer from the object domain 
into the image domain, but owing to limiting factors such 
as the finite size of the energy source, detector size, 

10 sampling frequency, display density, software filter 

function, and possibly partial-volume effects experienced 
with some imagers, an infinitely fine delta function in 
the object domain cannot be faithfully reproduced in the 
image domain. Instead, a smeared-out image, or point- 

15 spread function (PSF) , is observed. Similarly, an 

infinitely sharp edge-response function (ERF) in the 
object domain becomes a smeared-out ERF in, the image 
domain- When two image ERFs are close to each other, the 
smeared-out effects run into one another. The 

2 0 deterioration of the spatial resolution is often 

monitored by the percentage modulation transfer which is 
given by the ratio, expressed in percentage, of the 
amplitude of the modulation in the image domain to that 
of the object domain. 

25 

The smearing effect becomes more intense as the adjacent 
ERFs of discontinuities or contrast profiles get closer 

to each other (of as the spatial f requency of thfe 

modulation becomes higher) ; and this also causes the loss 

3 0 of profile heights. The modulation transfer function 

(MTF) is the function of the ^normalised' percentage 
modulation transfer over the whole range of spatial 
frequency until the cut-off frequency at which the 
amplitude of the image modulation transferred becomes 
35 zero, or, in practice, when it falls below the noise 
level and becomes undetectable. 
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The shape of the overall MTF over the full range of 
spatial frequency is similar to a roll-off sine function, 
and the high contrast profile has a higher cut-off 
spatial frequency than a low contrast profile. Within 
5 the working range before the cut-off spatial frequency, 

the inherent loss of the spatial resolution (that is, the 
part that is indicated by the smear ed-out effect on the 
corner edge of the ERF) cannot be restored or partially 
restored even by re-scanning the image with an ultra high 
10 resolution digital scanner system. 

It is an object of the present invention to provide a 
method and system by which the above problem can be 
overcome at least partially. 

15 

According to one aspect of the present invention there is 
provided a method of image enhancement wherein a de- 
convolution process is applied to the image-domain 
results of an object-scan to re-construct therefrom the 

20 respective point-spread function (PSF) effective in the 
object- to image-domain transfer at one or more object- 
discontinuities, and correlating this re-constructed 
function with the image-domain results of those one or 
more discontinuities for image-resolution enhancement as 

25 to discontinuity-location in the image domain. 



According to another of the aspects of the present 

invention there is provided an image-enhancement system 
comprising means for performing a de-convolution process 

3 0 on the image-domain results of an object-scan to re- 
construct therefrom the respective point-spread function 
(PSF) effective in the object- to image-domain transfer 
at one or more object-discontinuities, and means to 
correlate this re-constructed function with the image- 

3 5 domain results of those one or more discontinuities for 
/image-resolution enhancement as to discontinuity-location 
in the image domain. 



A least-squares running filter may be used for the 
de-convolution process. 

In the method and system of the present invention the 
5 optimum edge position corresponding to the discontinuity 
or sharp-edge position of the object ERF may be pin- 
pointed from the central point of the full-width half- 
maximum (FWHM) of the image PSF; the pin-pointing may be 
to sub-pixel accuracy for the image ERF. The original 

10 discontinuity or sharp-edge feature of the ERF may then 
be restored within the image domain by removing the 
sub-pixel values from outside the optimum edge position 
to compensate for those within. It is to be noted that 
the sub-pixels become pixels, and that the enhancement is 

15 equivalent to the performance of an extra high resolution 
image transfer system. 

It may be shown mathematically that the image ERF of an 
infinitely-sharp edge-response function ERF can be 

2 0 reproduced by convoluting the image PSF with the object 

ERF. This is illustrated in Figure 1 of the accompanying 
drawings. There is also illustrated in Figure 1 how the 
image ERF is used in conjunction with the image PSF to 
determine the optimum edge position from the full-width 

25 half maximum FWHM. The optimum edge position 

corresponding to the sharp edge position of the object 
ERF may be pin-pointed to sub-pixel accuracy for the 
image ERF from the central point of the FWHM. 



3 0 For one-dimensional cases, the operation in accordance 
with the invention is relatively simple, as only either 
the X or the Y profile (that is to say, as indicated 
below, a line spread function (LSF) ) is involved). But 
for the two-dimensional operations, both the X and Y 

35 profiles, and if necessary, the XY diagonal profiles to 
eliminate any possible streakings in the image, may be 
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involved. In this case, a proper weighting scheme will 
be required to re-construct the image, 

Once the original sharp-edge feature has been pin- 
5 pointed, that feature may be restored within the image 

domain by removing the sub-pixel values from outside the 
optimum edge position and utilising them for compensation 
of those within. This enhancement operation is 
illustrated in the form of the modulation display of 
10 Figure 2 of the accompanying drawings. 

The technique of the invention may partially recover the 
loss of the spatial resolution without the trade-off loss 
of other properties such as image noise. Furthermore, 

15 the enhancement of the spatial resolution in an image, of 
the invention may lead to enlargement or magnification of 
a small area or particular region of interest (ROI) 
without blurring (or step) effects at the profile edge, 
for example as when the sub-pixels are displayed as 

20 pixels at the ROI. 

A feasibility study using a theoretical data set and the 
noisy low-contrast magnetic resonance (MR) images of pig 
brain and lamb kidney surrounded by water have been 

25 successfully carried out. The results are illustrated in 
Figures 3 to 6 of the accompanying drawings, in which 

Figure 3 shows results for a theoretical data set, 

Figures 4 and 5 the results achieved for pig brain, and 
Figure 6 for lamb kidney. The signal-to-noise figures of 

3 0 the re-created PSF (measured by the main peak height to 
that of the adjacent peaks) were over 8/1. 
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Figure 7 of the accompanying drawings is a flow chart 
illustrative by way of example of the image-enhancement 
system. 



The method and system of the invention are applicable, 
for example, in medical science, in engineering, in 
physical science and in the field of instrumentation 
generally. In the context of medical science the 
5 invention is especially applicable to diagnostics and 
treatment planning, and may be utilised in connection 
with magnetic resonance (MR) imaging, computer assisted 
tomography (CT) , X-ray radiography, film or print copy or 
digital form, digital X-ray f luorography , ultra-sound 

10 imaging, nuclear medicine, positron emission tomography 
(PET) and other camera or imaging* In MR, X-ray, and CT 
it may be used, for example, to pin-point the true 
boundaries of tumour growth or any other abnormal effects 
essential for the diagnosis, surgery, and radiation 

15 treatment. Moreover, the technique is particularly 

suitable for use in X- ray digital f luorography , in which 
small structures under study are highlighted by injection 
of contrast liquids; the small structures may also be 
isolated from surrounding interfering effects by using an 

2 0 image subtraction technique. The present high-resolution 

enhancement technique may be used to improve the image 
quality further. 

Image subtraction coupled with the present resolution- 
25 enhancement technique may also be applied, for example, 
to MR images from two different scan-sequences. Also, 
the magnification of an image ROI without blurring (or 
step) effects at the edges of the profile has many uses. 

3 0 The inherent resolutions of X-ray radiography, 

ultra-sound imaging, nuclear medicine, and PET scanning 
are relatively low, and some are used for real-time 
study. Only the individual still frame or hard-copy 
images may be re-processed. 
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In the context of engineering, physical science and the 
field of instrumentation, the invention is applicable to 



one-dimensional imaging as used, for example, in regard 
•to bar-code patterns, the spectrum of DNA analysis, iris 
patterns of eyes (for example, for identification 
purposes in commercial banking) , finger-print 
5 identification, and emission spectroscopy. The invention 
is also applicable to two-dimensional imaging, for 
example, in relation to images obtained by satellite or 
pattern recognition, or from a surveillance camera or 
during laboratory experimentation. As a general matter, 

10 the invention is applicable where there needs to be 

accurate determination of the edge position in an image 
versus the true object-edge position, for the purpose, 
for example, of measurement of the positional 
displacement between object and image, distortion 

15 correction and manufacturing control. 

As an example of one of many applications of the present 
invention, a method and system that uses CT and MR 
imaging in conjunction with one another, will now be 

2 0 described. 

The major contribution to the magnetic-resonance (MR) 
signal comes from the abundant protons content of water 
molecules and protein. It is a quantum process at the 

25 Larmor frequency according to the magnetic field in use. 
The ^Tl-weighted' and ^T2 -weighted' MR signals from 

protons provide contrast numbers that are relative in 

scale, whereas in CT, the X-ray absorption is a 
polychromatic attenuation process affected by the 

3 0 electron densities of all the atoms presented within the 

X-ray beam. There is no equation to correlate the CT 
number (or the linear attenuation coefficient, electron 
density, or tissue density) with the MR-contrast numbers; 
no direct calibration between the two types of signal is 
35 possible. This lack of correlation is confirmed by 

consideration of bone and air which are at opposite ends 
of the CT contrast (absolute) scale using water as the 



base-line reference, but which are at the same end of the 
MR-image contrast (relative) scale owing to their conation 
low proton-population. 
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The lack of correlation between the CT and MR signals 
acts against their use in combination for imaging 
purposes, but the present invention provides a method and 
system by which the advantages of each may be utilised to 
improve image resolution and contrast information. 



In the latter regard, CT provides a high spatial 
resolution but only in regard to view normal to the 
transverse slice. Resolution for all re-constructed 
non-transverse planes is poor owing to the need to use 
15 elongate voxels to improve signal-to-noise ratio. Also, 
the partial-volume effect of using elongate voxels may 
give rise to detection errors at the thin edge of a 
contrast profile of a lesion. MR, on the other hand, can 
give the same high degree of spatial resolution viewed in 

2 0 the normal direction to any image-slice plane, and can 

also provide isotropic resolution with cubic-voxel volume 
imaging. 

To this end, multiple-slice transverse CT scans are 
25 collected across a section of the volume of interest in a 
patient or other subject. Corresponding multiple-slice 
transverse MR scans of the same volume are collected 
correspondingly. The slice thickness of the latter scans 
may be one half, or smaller, of the thickness of the CT 

3 0 slices, and may be collected two-dimensionally or three- 

dimensional ly. The patient is constrained throughout on 
a couch that provides landmarks with respect to a 
coordinate reference arrangement on the couch-top. This 
is to ensure the reproducibility of anatomical positions 
35 and features to the first order accuracy for the 

corresponding CT and MR scans, and also for the final 
radiation treatment to be made. 




The respective transverse planes of the CT and MR images 
are processed individually in the method and are matched 
with one another in a de-convoluted space for the CT and 
MR images, as later explained. The two sets of 
5 de-convoluted maps are then merged together to a second 
order of accuracy in order that the CT numbers may be 
transferred over to replace the corresponding MR contrast 
numbers. Once this has been achieved, non-transverse (or 
oblique) planes can be obtained from the two-dimensional 
10 (2D) MR images or from the re-arrangement of the 

corresponding voxels of the three-dimensional (3D) volume 
images; where a 2D-MR is used, a further step of contrast 
transformation may be required. 



15 In the method, CT multiple-slice transverse (Tr) scans 

are collected across a section of the volume of interest 
in a patient. Similarly, corresponding MR multiple-slice 
transverse scans of one half, or smaller, of the slice 
thickness of the CT slices are collected across the same 

2 0 volume of interest by 2D-MR sampling or 3D-MR sampling. 



The patient is constrained from movement during and 
between the two scanning processes on a couch-top (the 
^universal' couch-top described below with reference to 

25 Figures 14 and 15, or the ^alternative universal' couch- 
top described below with reference to Figures 24 and 25) , 

which also provides landmarks within a coordinate 

reference arrangement defined in the couch-top. The 
constraint and landmarking is necessary to ensure that 

3 0 anatomical positioning and features are maintained the 
same to a first order of accuracy for the CT and MR 
scans, for diagnostic purposes as well as for the final 
radiation treatment. The respective transverse planes of 
CT and MR images are processed individually by using 

35 ^boundary' or ^finger-print' matching techniques in a 
de-convolution space (DCS) for the CT & MR images. In 
these transverse image planes, in particular, the skin- 



contours features along the sides of the patient may be 
best used for second-order alignment and matching 
purposes, as they are less affected by patient -movements; 
see Figure 26. The transitional error, dx, may be 
5 readily corrected with respect to the coordinate 

positions of a rectangular tubing system embedded in the 
couch-top (for example as described below with reference 
to Figures 24 and 25). The processed data may then be 
used for a ^diagnostic and statistics software package^ 
10 of CT image versus MR image for their exactly 

corresponding transverse slice (s), and an associated 
^statistical package' for accurate computation of the 
^true' area, and then the 'true' volume, of a lesion or 
tissue profile or contour. 

15 

The two sets of de-convoluted maps may also be merged 
together to a second order of accuracy in order that the 
CT numbers may be transferred over to replace the 
corresponding MR contrast numbers. Once this has been 

2 0 done, non-transverse (or oblique) planes are obtained 
from the 2D-MR image or from the re-arrangement of the 
corresponding voxels of the 3D-MR volume; in the 2D-MR 
case, a further step of contrast transformation is 
required. The transferred contrast data may then be used 

25 in a ^3D-radiotherapy treatment planning software 

package' for an in-plane oblique-image pseudo-3D approach 
using MR images. 



Geometrical distortion in the MR images may be corrected 
3Q to a certain extent, by a combined spin-echo (SE) and 
field-echo (FE) scan sequence, with the SE part to 
correct for the field inhomogeneity in BO, and the FE 
part to overcome the impurity or non-linear gradients. 
The geometrical distortion may also be dealt with 
35 empirically using a post-process correction based on a 

set of phantom data relating to a drum phantom. Movement 
problems may be suppressed by some scan sequences and/ or 
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by gating; for certain parts of the body there is no 
movenient problem for the MR scan. 

Once the necessary sets of MR images with transferred CT 
5 numbers have been obtained, 3D-RTP may be started using 
the current 2D-RTP software package. The whole ^3D-' 
space is covered through the use of several 2D planes in 
an in-plane oblique-image pseudo-3D approach. 

10 The optimum-edge determination software sub-routine of 

the two-dimensional de-convolution technique (2D-DCS) may 
be used to help map out the ^true' boundaries of the 
various tissue profiles of an MR image. When the MR 
scans are carried out on the alternative universal couch- 

15 top described below with reference to Figures 24 and 25, 
embedded with the five MR solutions identified in Table I 
below, the set of Tl and T2 relaxation times and the 
standardized proton density with respect to a CUSO4.5H2O 
solution of Tl = T2 = 2 00 ms, may be each derived for the 
20 various tissue contours from the ratio of two MR signal 
intensities of two different MR scan sequences. They 
form the ^identification' set from which the tissue types 
may be identified. 

2 5 The results obtained may be used to produce an MR 

diagnostic and statistics software package. From the 

identity of the tissue type and through the Bulk 

Heterogeneity Correction method (described by Richard A, 
Geise et al Radiology, 124:133-141, July, 1977), the 

3 0 normalised electron density may then be assigned to each 

tissue type. With all the tissue contours accordingly 
mapped out with their respective electron densities, the 
radiotherapy treatment planning may be smoothly carried 
out on any MR image. The actual oblique MR image 
35 planning will lead to the eventual software package of 
^3D-radiotherapy treatment planning' using MR images 
alone. Before the actual planning is started, any 
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remaining geometry distortion within an MR image may be 
empirically and accurately corrected by a DCS sub-routine 
using the geometry distortion correction data in-stored. 

5 The matching technique is based on a method using de- 
convolution space (DCS) in * combination with a least- 
squares curve-fitting method. These methods are 
separately well known in pure mathematics, but not, as in 
the present method and system, combined to produce de- 
10 convolution maps with the least-squares fitting curve 

used as a running filter. The filter is used to overcome 
practical difficulties due to noise within the image 
pixels. 

15 For a mathematical model of the method, consideration is 
given to a line spread function (LSF) arising in the 
imaging system from an impulse or S function; this is 
represented in Figure 8 of the accompanying drawings. 

2 0 An example of the way in which the LSF of Figure 8 is 

convolved with an object step-down function to result in 
a roll-down edge response function (ERF) , and how this 
through the de-convolution process reproduces the LSF, is 
illustrated in Figure 9 of the accompanying drawings. As 

2 5 illustrated, the LSF is a one dimensional profile across 
the edge of a contour with its peak position representing 
the 50%-level point of the ERF. 



Figure 9 and the mathematical expressions set out below, 
indicate the way in which the LSF can be recovered from 
the resulting ERF at a contrast edge. For convolution 
related to points 1 to n of the step function: 

at point 1: ^ 

J^LCx) Ax = E(x)^ 



35 
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at point 2: 



a-l 



at point 3 : ^ ^ 

Ax - f^L(x) Ax = £;(x) 3 



a-2 



5 until at point nj 

J^Lix) Ax - J^L(x) Ax = E(x)^ 



For de-convolution ie, for the ERF to manifest the LSF: 

a 

E(x)^ - E{x)^ - 5^L(x)Ax = L(x) (3.^)^3 Ax 



a-l 



{^(x)i - £'(x)2)/Ax = i:(x) (3.1^^3 



10 thus: 



{^(x), - £'(x)3)/Ax = L(x) 



(a-2)-{a-l) 



leading to 



Thus, the same shape of LSF is recovered by the 
15 de-convolution process independently of the sense of the 
ERF; it does not matter whether a roll-up or a roll-down 
ERF of a positive or a negative contrast profile, is 
involved. This is an important property as a positive 
contrast contour will have roll-up (from low to high) 
2 0 ERFs at both ends whereas a negative contrast contour 

will have roll-down (from high to low) ERFs at both ends. 
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In the practical implementation, the LSF is derived from 
the ERF by de-convolution using a running filter. This 
enhances the accuracy of the method in overcoming the 
problem of noise that affects the digitised pixel values 
of the image. Use is made of a least-squares fitting of 
a set of data points that advances along the whole length 
of the function from one point to another. 



Assuming that y == aQ + ajX + ajX^ represents the ERF curve 
10 a five-point or seven-point fit is used, and the normal 
equation becomes: 



El Ex^ Ex| 
Jlxj^ Hxf Exi 
Exj Ex"j Ex^ 



aO 
al 
a2 



where: i is 1+n, ... 5+n; ... until m-5, ... m 
or 

15 i is 1+n, ... 7+n; ... until m-7, ... m 

n is 0,1, 2, ... ; and 

1, ... m is the span of the ERF profile. 



The solution may be derived from either: 



aO 


El 


Ex 


Ex2 


-1 


Ey 


al -. 


= Eat 


Ex2 


Ex^ 




Hxy 


a2 


Ex2 


Ex3 


Ex* 
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or : 



aO 



Ey Ex Ex 2 
Exy Ex^ Ex^ 
Ex2y Ex3 Ex* 



al 



El Ey Ex 
Ex Exy Ex' 
Ex^ Ex^y Ex* 



a2 



El Ex Ey 
Ex Ex^ Exy 
Ex^ Ex' Ex^y 



El Ex Ex2 
Ex Ex2 Ex' 
|Ex2 Ex' Ex* 



For both of these equations to be valid:- 
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El Ex Ea:^ 
Ex Ex2 Ex3 
|Ex2 Ex^ Ex^ 



^ 0 



The gradient at dy/dx^^.„) or dy/dX(4+^) can then be derived 
and plotted against x for the LSF profile. 

5 The graph of dy/dx against x gives the LSF profile. The 
peak of this profile is located half way between the mid- 
points of the ascending and descending limbs of the graph 
and is derived from this as the so-called full-width- 
half -maximum (FWHM) with an accuracy to the sub-pixel 
10 level. 

The point spread function (PSF) is the two dimensional 
profile which may be derived, in practice, from the two 
corresponding LSFs orthogonal to one another within an 
15 image plane. The peak position of the PSF profile is, 
therefore, from the ^mean' or ^cross-over' of the two 
peaks or the two LSF profiles. The PSF is obtained in 
practice from two orthogonal axes in a two-dimens;ional 
plane. 

20 

The waveform relationships between real space (RS) , image 
space (IS) , and de-convolution space (DCS) are 

illustrated at fa) , fb) and (c) respectivelv of Figure 10 

of the accompanying drawings. These show the 

25 relationship between peak positions pi, p2 , pi', p2',pl" 
and p2" and dip positions d' and d" in the DCS, with the 
positions and magnitudes of contrast value of the 
waveforms in the image space (IS) and real space (RS) , 
and highlight certain useful properties of the process. 

30 In particular, the tail-end to tail-end distance of each 
LSF profile gives a spread of 100% down to 0%, and vice 
versa, in the corresponding ERF profile in IS. Also, the 
half way points of the full-width-half -maximum (FWHM) in 
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the LSF profiles give the peak positions of pi, p2 , pi', 
p2 ' , pi", and p2" which correspond to the 50% level 
points of their related ERF profiles; the dip positions 
d' and d" in the LSF profiles correspond to the central 
5 maxima or minima of the profiles in IS, and the distances 
of pi' to p2 ' , and pi" to p2" correspond to the FWHMs of 
the profiles in the IS and the distances between the 
edges of the step function in RS. 

10 From the above, the generation of the LSF (or PSF) is 

independent of the roll-up or roll-down nature of ERFs at 
the edges of the contrast contour. In other words, it is 
independent of the sense and the absolute value of the 
contrast numbers within the contour. The peak position 

15 of the LSF (or PSF) is the central half-way (50%) point 
of the roll-up or roll-down ERF which is the optimum 
position at which the two sets of ERFs are matched 
together with the overall effect of not distorting the 
contour shape within the boundary of the ERFs. 

20 

In order to match or merge together the two corresponding 
CT and MR images, they are transformed into DCS. In DCS 
the peaks of the LSFs (or PSFs) mark the central half-way 
points of the ERFs of the original two-dimensional image 

25 planes. These peak positions are then linked together to 
form the Me-convolution' or ^finger-print' maps. 
Through the assistance of these maps, the two image files 
may now be matched or merged together to a second order 
of accuracy. Figures 11 and 12 of the accompanying 

3 0 drawings illustrate one-dimensional and two-dimensional 

generation of de-convolution or finger-print maps in DCS. 



The subsequent transformation of the contrast numbers is 
guided by a moving sub-set filter conformed by the 
3 5 sub-set mean values and the overall uniformity of the 

profile of the MR image; this filter is quite similar to 
the running filter used for the conversion of ERF to LSF. 



• 
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The number of samples in the moving sub-set filter may 
also be used for the cut-off rejection of small or 
^ point' contours, that is to say, spurious contours 
introduced by noise, within the CT or MR image. 

5 

A similar technique is employed for MR(Tr) to MR (Oblique) 
transformation of contrast numbers. Since a high spatial 
resolution and small partial volume effect are 
maintained for the thin 2 D-MR (Oblique) slice or the 
10 3 D-MR (Oblique) plane, the transformation is relatively 
simple. Through this DCS technique, the high spatial 
resolution operation is maintained throughout the whole 
process of re-processing . 

15 The combined MR and CT method of the invention requires 
the matching of corresponding CT and MR images with one 
another, and it is important that the same ^laid-out' 
features of anatomy are presented for the each scanning 
activity, and that the same features are presented again 

2 0 for radiation treatment. The couch-top is used for this 

and the landmarks and reference axes defined in it are 
encoded into the images, for use as guides for the 
repositioning of the patient for any subsequent scanning 
or other activities, and also for the first order 

25 alignment in the initial matching of CT and MR images 
(allowing more refined matching to the second order of 

accuracy using the DCS technique, to be initiated) . 

The form of the couch-top used will now be described with 

3 0 reference to Figures 13 to 15 of the accompanying 

drawings. 

Referring to Figure 13 to 15, the couch-top has a flat 
upper part 100 that is carried by a curved lower part 101 
3 5 which fits to the scanner table of the CT or MR 

installation. The patient lies on the flat part 100 and 
is constrained on it (by an arrangement not shown) , so as 
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to reproduce the same anatomical features for scanning 
and treatment; the treatment table similarly has a flat 
surface. 



An X-tubing array 102 is embedded into the upper part 
100. The tubing of the array 102 is of plastics or glass 
having an external diameter of up to 6 mm and an internal 
diameter of up to 3 mm, and is filled with a chemical 
solution of MnCl2. The cross-sectional image of such 
tubing will show a positive profile for MR and a negative 
profile for CT. However, the DCS technique of processing 
these images pin-points their centres irrespective of the 
contrast values represented. 

15 As shown in Figure 14, the patient is constrained by the 
constraint arrangement to lie lengthwise of the flat part 
100 (ie along the z-axis of the coordinate system defined 
in the couch-top) and centrally of it (ie with respect to 
the lateral x-axis) ; the patient's head is placed on a 

2 0 head-rest 103 suitable for the scan, CT or MR, to be 

undertaken. Three landmarkers/ constraint levers 104, two 
under the armpits and one at the crotch between the 
thighs, are also used to position the patient as well as 
to provide some form of conformation with respect to the 




2 5 reference coordinate system. Air-compressed cushions 



(not shown) may be gently pressed in from either sides to 
conform the limbs and hence the body along its length. 




The X-tubing array 102 inside the bottom flat part 100 
3 0 may be used as the reference system for the transverse 

and sagittal planes, and the oblique planes between them. 
The flat surface of the part 100 or the flat contour of 
the patient's skin surface in contact with it may be used 
as the datum for the coronal planes. However, an extra 
3 5 small X-tubing array (not shown) may be placed above the 
patient using the three landmarkers as the supporting 
legs; this enables the small-angled oblique planes 



18 



between the coronal and transverse planes to be 
monitored. 

If the internal-bore diameter of the magnet of the MR 
5 scanner is large enough, an extra X-tubing array 105 

(Figure 15) may be used along either side of the patient 
to help in monitoring the small-angled oblique planes 
between the coronal and sagittal planes. Furthermore, a 
curved X-tubing array (not shown) may be placed inside 
10 and on the bottom of the head rest (or head coil) 103 
used for the MR scan for the purpose of clearly 
establishing the coordinates of the head. 

If the internal-bore diameter of the magnet of the MR 
15 scanner is not large enough for the extra array 105 to be 
fitted, then the ^alternative universal couch-top' of 
Figures 24 and 2 5 may be used; Figure 2 5 is a section 
taken on the line XXV-XXV of Figure 24. As illustrated 
in Figures 2 4 and 25, the bottom X-tubing pattern is in 
20 this case replaced by tubing 110 having a zig-zag pattern 
throughout the length of the couch-top. Extra capacity 
is provided by straight tubing 111 that extends centrally 
lengthwise of the couch-top. 

2 5 The tubing 110 of zig-zag pattern may be as illustrated 

in the detail of Figure 24 and in Figure 25, formed by 
double-layer tubing of rectangular cross-section, so that 
edge-position determinations may become one-dimensional 
when applying the de-convolution space subroutine (i.e. 

3 0 ID-DCS subroutine) ; the tubing 111 may of single-layer 

rectangular cross-section as illustrated. The 
rectangular tubing 110 and 111 may be filled as indicated 
in Figure 2 5 with the five MR solutions to S4 of Table 
I, for standardisation and calibration purposes for 
3 5 applications in the MR ^diagnostic and statistics 

software package' and the stand-alone MR method for the 
^3D-radiotherapy treatment planning software package'. 
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From the five MR solutions Sq to S4, choice may be made 
from the four MnCl2.4H20 solutions, to S4, to cover the 
full range of values of Tl and T2 for anatomical tissues, 
and of the fifth solution, Sq, of CUSO4.5H2O, that is 
^equivalent' to "loosely bound water" (see pages 34-38 of 
^NMR, A Primer for Medical Imaging' by G.L. Wolf et al , 
published by Slack Incorporated, 1984) • 



10 



Examples of solutions Sq to S4 are as set out in Table I, 
as follows: 

Table I 



Solution 




at 0.5T 


T2 

at 0.5T 


So 


1.2 5 g/1 CUSO4.5H2O 


200 ms 


200 ms 


Si 


3.41x10^** Mn**'^ ions/ml 


84 0 ms 


3 00 ms 


S2 


1.15x10^'' Mn"*"^ ions/ml 


44 0 ms 


120 ms 


S3 


2.30x10^'' Mn+2 ions/ml 


2 50 ms 


60 ms 


S4 


4.37x10^'' Mn"*"^ ions/ml 


150 ms 


3 0 ms 
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Figure 27 shows a form of cover that may be used for 
suppressing movement of, in particular, the patient's 
abdomen. The cover is of polystyrene sheet 12 0 with 
cushioning 121 around apertures for the patient's legs. 
The sheet 12 0 is secured over the patient and between 
25 his/her legs using fabric elements 122 that establish 

with an underlay (not shown) releasable fixings of fabric 
hook-and-eye form (for example, as sold under the 
Registered Trade Mark VELCRO) • 

3 0 The material of the couch-top is carbon fibre so as to 
ensure a low linear attenuation coefficient with 
transparency to X-rays. Edges of the couch-top are 
smoothly rounded to avoid sharp edges that would cause 




20 

streaking in the CT image. Non-magnetic material is used 
throughout to avoid disturbing the uniform magnetic 
field, BO, used for MR. 

5 The radiation treatment table needs to be capable of both 
trans lational, "in and out" movement as well as 
rotational movement to cover the 3D planning and 
treatment. 

10 A final empirical correction of any residual geometrical 
distortion may be carried out using a cylindrical or 
elliptical drum phantom. A cylindrical drum phantom is 
illustrated in Figures 16 and 17, and an elliptical drum 
phantom in Figures 18 and 19 of the accompanying 

15 drawings. 



Referring to Figures 16 and 17, the cylindrical drum 
phantom, which may be of about 160 mm in length and is 
for use with the head coil of the scanner, includes a 
central rod 200 of polymethyl methacrylate (as sold under 
the Registered Trade Mark PERSPEX) having a diameter of 
about 2 mm and defining the phantom centre, Cph. The 
body of the phantom is built up as concentric sections 
about the rod 2 00, each comprising evenly-spaced blocks 
2 01 of polymethyl methacrylate and with the blocks 2 01 of 
each section in exact register radially with the gaps 202 
between the blocks 201 of the preceding section. The 
central region 203, having a diameter of about 20 mm, may 
be left blank (since the magnetic field within the 
30 central region is uniform to within 0.1 ppm, so no BO 
correction is required) . The gaps 202 and the central 
region 203 are filled with an NMR solution (e.g. MnClj or 
CUSO4) . 

35 Peak positions of the PSFs in the DCS map for the drum 

are derived and stored. This data is then corrected for 
each point of the map with respect to the iso-centre of 
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25 



10 
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the scanner system, Ciso, using the following vector 
equation:- 

Ar = rj - r, 

where Ar is the correction vector between Ciso and 

Cph ; and 

r2 and ri are the vectors to the relevant point 
of the map from Ciso and Cph respectively. 



The elliptical drum-phantom of Figures 18 and 19 is 
constructed in essentially the same manner as the 
cylindrical drum phantom, and is for use with the body 
coil of the scanner. The major and minor axes are chosen 
15 according to the field of view being studied. 

As an alternative,, the cylindrical and elliptical drum 
phantoms may be constructed as shown in Figures 2 0 and 
21, using thin plates 300 of polymethyl methacrylate 

2 0 instead of the blocks used for the phantoms of Figures 16 

to 19. The plates 3 00, which may have a thickness of 
between 1.5 mm and 2.0 mm, are spaced radially from one 
another by gaps 3 01 that, together with the central 
region 3 02 are filled with appropriate NMR solution. 

25 

CT scan images are not subject to geometrical distortion 
so CT scan data accumulated using the phantoms, and the 
DCS profile data consequently derived and stored, are 
used as a datum or reference set for correction of the MR 

3 0 scans. Such corrections also automatically correct for 

manufacturing tolerances, and the DCS data derived from 
the phantom through the CT and MR scans represents a 
reference situation devoid of motion. 



35 



The DCS technique described above is applicable further 
in diagnosis and treatment. In particular, it may be 
used to assist visual detection of gradual change of 
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contrast level in scan images. For example, as 
illustrated in Figure 22 of the accompanying drawings, 
the 100% point at which the roll-off of a lesion profile 
begins can be readily plotted from the DCS profile, 
5 together with the 50%-level point, and then the 0% point 
where the profile has completely levelled off- Graphic 
representations of these various points can then be 
plotted out to facilitate diagnosis, as illustrated in 
Figure 23 of the accompanying drawings, in the form of a 
10 string of characters 400 delineating the 100%-level, a 

string of characters 401 delineating the 50%-level and a 
string of characters 402 delineating the 0% level. 

In current radiation treatment planning, a composite 
15 volume leaving a margin of 15 mm around a lesion, is 

used. To be effective, the tail-end of the lesion growth 
is to be within the treatment-beam path, and in this 
regard there is advantage in extending the DCS-profile 
plot of Figure 2 2 to include a point which is a distance 
20 (d) beyond the 0%-level. This is then used for diagnosis 
and treatment, to establish an outer string of characters 
403 in the graphical representation of Figure 23. The 
distance (d) covers the penumbra-spread as well as any 
possible mis-alignment by the operator. 

25 

The de-convolution space (DCS) technique leads to two 

major software packages for: (a) optimum edge 

determinations, and (b) spatial resolution enhancement. 
Applications that may be derived from these packages are 
30 identified in the flow-chart of Figure 28. 



From the software for both optimum edge determinations 
(for image profiles and artificial edges) and spatial 
resolution enhancement for all one- and two-dimensional 
35 images, products for both medical and non-medical 
applications may be provided as follows: 
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(1) software package for the spatial 
resolution enhancement and the enlargement of 
all two-dimensional image-regions of interest 
free of smearing (or step) effects at the 

5 image-profile edges; 

(2) software package for the spatial 
resolution enhancement and enlargement of all 
one-dimensional images-regions of interest free 

10 of smearing (or step) effects at the image- 

profile edges; 

(3) software package for accurate edge 
determinations for all one-dimensional and 

15 two-dimensional images for multiple 

applications; and 

(4) image-enhancement hardware system to cover 
the areas of software packages (l)-(3). 

20 

From the software for optimum edge determinations (for 
image profiles and artificial edges) alone, the following 
hardware and software packages for medical applications 
may be provided: 

25 

(1) a section or full-scale section of the 
^alternative universal couch-top' design for 
first- and second-order alignments and 
standardisation purposes for applications as 

3 0 identified in items (2) to (6) below; 

(2) an empirical geometry-distortion 
correction package for MR imaging; 



35 



(3) a diagnostic and statistics package of 
2D(CT) image versus 2D(MRI) image of exactly- 
corresponding transverse slices - the 
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Statistics package may help to compute the 
accurate area within the ^true' boundaries of a 
lesion or a tissue profile, and then follow-up 
with the eventual calculation of the ^true' 
5 volume of such profile or effect (other 

statistical approaches may also be performed 
within the tissue boundaries - for example, 
within these boundaries, noise grains may be 
ironed out if necessary by taking the mean of 
10 the output of a long-running filter across the 

whole profile) ; 

(4) MR diagnostic and statistics package to 
compute accurately the Tl and T2 relaxation 
15 times and the standardized proton density of 

the tissue profiles for diagnostic purposes, 
and to perform the same statistics as that of 
package (3) ; 



20 
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(5) a combined CT & MR method using the 
multiple-slice data of item (3) for 
three-dimensional radiotherapy treatment 
planning within the oblique MR image planes; 
and 

(6) a stand-alone MR method using the 
calculated Tl, T2 and the standardised proton 



density within the tissue profiles of the 
multiple-slice image-data of item (4) in order 

3 0 to identify the tissue types and then to obtain 

the electron densities through the Bulk 
Heterogeneity Correction method so that the 
three-dimensional radiotherapy treatment 
planning within the general oblique MR image 

35 planes may be carried out. 
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Fig. 22 Fig. 23 
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